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ABSTRACT 


As  timber  design  codes  move  from  allowable  stress 
design  methods  to  limit  states  design,  there  is  need  to 
define  the  behaviour  of  timber  members  up  to  ultimate 
conditions.  The  present  study  is  an  investigation  into  the 
behaviour  and  ultimate  strength  of  commercial  grade 
Douglas-fir  glued  laminated  timber  beam-columns. 

Approximate  analytical  procedures,  based  on  an 
elasto-plast ic  compression  and  linear  tension  stress-strain 
distribution  in  tension,  are  used  to  predict  ultimate 
strength.  Interaction  curves  based  on  material  properties 
obtained  from  small-scale  tests  are  developed  from  the 
analyses . 

The  experimental  program  consisted  of  testing  nine 
full-scale  factory  manufactured  beam-columns,  ten 
compression  specimens  and  ten  standard  small-scale  tension 
specimens.  Variables  investigated  in  the  beam-column  tests 
included  slenderness  ratio  and  magnitude  of  axial  load.  The 
behaviour  of  the  beam-column  specimens  was  monitored  by 
measurements  of  lateral  loads,  cross-sec t ion  strain 
distribution  and  lateral  deflections. 

Interaction  curves  derived  from  Newmark's  numerical 
analysis  procedure  and  modiified  by  an  undercapacity  factor 
of  0.7  are  in  good  agreement  with  test  results. 
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1.  INTRODUCTION 


1 . 1  General  Remarks 

Glued-laminated  timber  columns  are  frequently  subjected 
to  situations  where  continuity  conditions  and  effects  of 
lateral  forces  such  as  wind  impose  significant  moment  in 
addition  to  axial  loads.  Presently  in  Canada,  timber 
beam-columns  are  designed  according  to  the  Code  for 
Engineering  Design  in  Wood,  CAN3-086-M80 1 ,  which  is  an 
allowable  stress  code.  The  current  move  towards  a  limit 
states  code  requires  an  understanding  of  the  behaviour  of 
these  members  up  to  ultimate  conditions. 

A  number  of  investigators  have  focussed  on  the  strength 
of  timber  columns  subjected  to  axial  load  with  small 
eccentricities.  Relatively  few  studies  have  related  to  the 
interaction  diagram  approach  to  design  of  beam-columns  based 
on  ultimate  conditions. 

1.2  Object  and  Scope 

The  main  objectives  of  this  investigation  are: 

1.  To  develop  analytical  procedures  which  predict  ultimate 
strength  of  timber  beam-columns. 

2.  To  carry  out  preliminary  testing  on  commercial  grade 
glued-laminated  timber  beam-columns  to  observe  their 
behaviour  and  ultimate  capacity. 

3.  To  establish  a  basis  for  further  studies  on  ultimate 
strength  of  timber  beam-columns. 
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2.  LITERATURE  REVIEW 


2.1  Previous  Research 


Comparatively 

little 

research 

has 

been 

conducted 

timber  members  subjected  to 

combined 

axial 

load 

and  bend 

moment.  In  1954, 

Pearson 

2  invest 

iga ted 

the 

ef  f ects 

species,  slenderness  ratio,  eccentricity  of  load  and 
orientation  of  growth  rings  on  the  strength  of  solid  timber 
columns.  Over  400  specimens  were  tested,  including  more  than 
250  eccentrically  loaded  columns.  Slenderness  ratio  (L/d) 
ranged  from  5  to  50  with  eccentricity  to  depth  ratio  of  0  to 
0.5.  Based  on  the  test  results,  a  modification  of  Jezek's 
formula  (through  Pearson,  1954)  was  proposed  for  calculating 
the  maximum  strength  of  eccentrically  loaded  columns. 
Pearson  also  compared  his  test  results  with  the  predictions 
based  on  the  modified  secant  formula.  The  test  results  were 
found  to  be  in  much  closer  agreement  with  the  modified 
secant  formula  when  the  eccentricities  were  small  than  when 
they  were  large. 

Wood3,  in  1961,  presented  formulae  for  calculating  the 
safe  capacity  of  timber  columns  with  lateral  loads  and 
eccentric  axial  load.  The  formulae  were  essentially  those 
developed  by  Newlin  (through  Wood,  1961),  but  included  some 
extensions  and  applications  to  design  problems.  Wood  assumed 
that  the  maximum  stress  developed  under  combined  load  was 
equal  to  the  flexural  strength.  For  columns  having  a 
slenderness  ratio,  L/d,  of  11  or  less  (short  members),  an 
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interaction  equation  for  the  general  case  of  eccentric  axial 
load  and  lateral  load  was  derived.  For  long  columns,  having 
L/d  of  20  or  more,  a  model  column  with  initial  cosine  wave 
curvature  was  used.  By  assuming  a  small  curvature,  a  general 
interaction  equation  for  this  case  was  also  derived.  It  was 
then  suggested  that  a  linear  interpolation  be  used  for 
calculating  the  capacities  of  columns  in  the  intermediate 
slenderness  range  i.e.  L/d  between  11  and  20. 

In  1970,  Hammond,  Curtis,  Sidebottom  and  Benjamin5 
studied  the  effect  of  eccentricity  and  end  restraints  on  the 
strength  of  timber  columns.  The  stress-strain  diagram 
obtained  from  simple  tension  and  compression  as  shown  in 
Figure  2.1(a)  was  idealized  as  elasto-plast ic  in  compression 
and  linear  in  tension  as  in  Figure  2.1(b).  Expressions  were 
derived  for  the  axial  load  P,  and  bending  moment  M, 
expressed  as  functions  of  the  cross-section  and  the  depth  of 
yielding.  Interaction  curves  were  then  obtained  for  various 
depths  of  yielding.  Curvature,  at  a  given  section  was 
determined  from  the  strain  distribution  in  terms  of  bending 
moment  and  depth  of  yielding.  By  assuming  the  deflected 
shape  as  a  sine  curve,  and  linear  elastic  end  restraints, 
the  expression  for  the  eccentricity  was  obtained.  The 
collapse  load  was  then  obtained  by  trial  and  error,  given 
the  eccentricity  and  the  depth  of  yielding  as  shown  in  Fig. 
2.2.  Theoretical  curves  relating  end  constraint,  column 
eccentricity  to  depth  ratio,  column  length  to  depth  ratio, 
and  average  fibre  stress  were  derived  from  the  analysis. 
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thus  simplifying  the  design  of  such  columns.  Hammond  et  al 
tested  54  solid  timber  columns  with  eccentricity  to  depth 
ratio  of  0.1  to  10,  three  different  elastic  end  restraints, 
and  slenderness  ratio  of  about  80  to  200.  The  test  results 
were  reported  to  be  in  good  agreement  with  the  theoretical 
curves . 

Zakic9,  in  1975,  provided  mathematical  solutions  for 
timber  members  subjected  to  bending  moment  plus  compressive 
axial  force  and  bending  moment  plus  tensile  axial  force.  He 
derived  interaction  equations  for  elastic  and  inelastic 
behaviour  in  both  cases.  Parabolic  compression  and  linear 
tension  stress-strain  curves  shown  in  Figure  2.3,  derived  in 
his  earlier  research6,  were  used.  Zakic  made  use  of  certain 
limiting  stresses  at  impending  yield,  i.e.  ayc=0.5auc  where 
ayc  is  the  compression  stress  at  yield  and  au c  is  the 
ultimate  compression  stress.  Adopting  the  mathematical 
equations  for  the  stress  diagram  proposed  by  Moe  (through 
Zakic,  1979),  non-linear  interaction  equations  were 
developed  for  inelastic  beam-column  behaviour.  Zakic 
reported  15  full-scale  glue-laminated  beam-column  tests.  The 
test  specimens  were  simply  supported,  and  two  symmetrically 
placed  concentrated  lateral  loads  were  applied.  The 
interaction  curve  obtained  from  the  test  results  was  in  good 
agreement  with  the  theoretically  predicted  curve. 

In  1979,  Larsen  and  Theiglaard10  derived  theoretical 
capacities  for  the  general  case  of  beam-columns  subjected  to 
equal  end  moments,  with  initial  double  curvature  and  initial 


' 


. 


' 

£  '  i  'l 

. 


1 1 


1 1  <j  i  tt s  i  <"  «  t 

■ 


■  * 


5 


torsion.  The  initial  curvature  and  torsion  were  assumed  to 
be  represented  by  cosine  functions;  and  the  failure 
condition  was  assumed  to  be  a  linear  interaction  of 
compression  and  bending  stresses.  These  assumptions  were 
used  to  solve  the  general  differential  equations,  leading  to 
an  expression  for  the  total  loading  capacity.  The  general 
expression  was  later  specialized  for  cases  of  combined  axial 
load  and  in-plane  bending  moment  or  pure  in-plane  bending 
alone.  Larsen  and  Theiglaard's  experimental  program 
consisted  of  tests  of  39  specimens  employing  different 
material  grades,  cross-sections  and  lengths.  For  the  case  of 
in-plane  bending  and  axial  load,  plumb-line  measurements  of 
initial  displacement  and  stress  ratios  from  the  Danish  Code 
(through  Larsen  et  al,  1979)  were  used  to  obtain  the 
ultimate  moment.  Close  agreement  between  the  measured  and 
calculated  ultimate  moments  was  reported.  Approximate 
expressions  and  design  curves  were  presented  for  use  in 
designing  members  for  the  cases  investigated.  It  was  noted 
that  column  capacity  was  not  only  affected  by  the 
slenderness  ratio,  but  also  by  the  depth  to  breadth  ratio; 
and  that  the  capacity  was  least  for  square  columns. 

Malhotra’1,  in  1982,  developed  a  mathematical  model  for 
the  analysis  of  timber  members  subjected  to  compression 
loads  with  small  eccentricities.  The  analysis  was  based  on 
Jezek's  simplification  of  an  ideal  elasto-plast ic 
stress-strain  relationship  in  compression  and  a  linear 
relationship  in  tension.  The  axis  of  the  deflected  column 
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was  assumed  to  take  the  form  of  a  half  sine  wave  . 
Equilibrium  conditions  were  established  at  the  column 
mid-height  to  obtain  the  critical  column  stress  as  a 
function  of  slenderness  ratio  and  eccentricity  to  depth 
ratio.  Malhotra’s  test  program  involved  over  350  columns 
with  eccentricity  to  depth  ratio  from  0.1  to  0.3,  and 
slenderness  ratios  ranging  from  40  to  160.  Most  columns  were 
pin-ended,  but  a  few  were  flat-ended.  The  test  results  were 
compared  with  the  modified  secant  equation  and  the 
Per ry-Robert son ' s  formula.  It  was  found  that  the 
experimental  results  were  in  good  agreement  with  Malhotra's 
model  based  on  Jezek's  approach  and  also  with  the  modified 
secant  formula.  Malhotra  proposed  a  semi -rat i onal  approach 
for  the  design  of  columns  subjected  to  axial  load  and  small 
eccentr ic i t ies . 

2.2  Present  Code  Requirements  and  General  Comments 

CAN3-086-M80  Code  for  Engineering  Design  in  Wood1  is  a 
working  stress  design  code.  Its  requirements  for  members 
subjected  to  combined  axial  load  and  bending  moment  are 
based  on  a  linear  interaction  equation  given  by: 

P/A  +  M/S  <1  (2.1) 
a  F 

where  P  is  the  concentrated  axial  load,  A  is  the  area  of 
cross-section,  a  is  the  allowable  unit  stress  in  tension  or 
compression  that  would  be  permitted  if  axial  load  only 
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existed,  taking  into  account  the  slenderness  ratio  and  the 
duration  of  load  factor.  M  is  the  total  bending  moment, 
including  the  moment  due  to  axial  load,  S  is  the  section 
modulus ;  F  i s  the  allowable  unit  stress  in  bending  that 
would  be  permitted  if  bending  only  existed,  and  modified  for 
the  duration  of  load. 

Some  researchers 3  9  -  1 1  -  ,  have  suggested  that  a  linear 
interaction  equation,  based  on  the  ultimate  strength  of 
timber  beam-columns,  is  conservative.  Larsen  and 
Theiglaard  10  reported  slightly  unconservative  results  when 
comparing  the  linear  interaction  diagram  with  strength  based 
on  assumed  initial  displacements.  However,  considering  the 
variation  in  the  properties  of  structural  timber,  the 
interaction  was  considered  acceptable. 
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F i gure  2 . 1 


Stress-Strain  Diagrams 
Compression 


in  Direct  Tension  and 


(a)  Cross  Section  (b)  Strain  Distribution  (c)  Stress  Distribution 


Figure  2.2  Strain  and  Stress  Distributions 
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3.  ANALYSIS 


3 . 1  Introduction 

Three  analytical  procedures  may  be  adapted  to  the  case 
of  timber  beam-columns  of  rectangular  cross-section;  namely 
the  method  of  assumed  deflected  shape  of  beam-column, 
Newmark's  numerical  integration  procedure,  and  the  moment 
magnifier  method.  The  following  assumptions  are  common  to 
the  the  first  two  methods: 

1.  Strains  are  linearly  distributed  across  a  section 
subjected  to  bending  moment  and  compressive  axial  load. 

2.  The  stress-strain  characteristics  for  wood  in  direct 
compression  and  tension  are  available  for  a  given 
moisture  content. 

3.  The  stress-strain  diagram  in  compression  is  idealised  as 
elat ic-per f ectly  plastic. 

4.  Modulus  of  elasticity  is  the  same  in  tension  and 
compression . 

5.  Failure  criterion  is  based  on  attainment  of  the  ultimate 
strain  in  compression  at  the  extreme  compression  fibre 
of  the  section.  Thus  the  failure  sequence  is  plastic 
compression-splitting  tension. 

6.  The  behaviour  of  the  beam-column  is  not  affected  by 
strength  reducing  defects. 

7.  Axial  load  and  bending  moments  are  in  the  plane  of 
bending . 

8.  Shear  forces  are  neglected. 
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9.  Columns  are  initially  straight. 

Figures  2.1  shows  typical  stress-strain  diagrams  in 
tension  and  compression,  and  the  idealization  of  the 
compression  stress-strain  diagram.  The  stress-strain 
distribution,  and  the  nomenclature  used  are  shown  in  Figure 
2.2. 

Based  on  the  above  assumptions  and  a  suitable  deflected 
shape,  interaction  diagrams  for  any  depth  of  yielding, 
slenderness  ratio,  and  ultimate  strains  in  compression  and 
tension  can  be  generated. 

3.2  Method  of  Assumed  Deflected  Shape  (Method  1) 

To  derive  the  interaction  equation  by  this  method5,  it 
is  further  assumed  that  the  beam-column  deflects  in  the  form 
of  a  portion  of  a  sine  curve.  Referring  to  Figure  2.2,  if 
the  axial  thrust  acting  on  the  member  is  P,  then 


P  =  J”A  adA 

or 

P  =  b[ayc  +  1/2  <7^  ( d-a-a  '  )  ] 

=  bdayc  +  0 . 5b ( d-a ) ( o yc +o u t ) 

where  a yc  is  the  yield  stress  in  compression,  out  is  the 
ultimate  tensile  stress,  a'  is  the  depth  of  yielding,  a  is 
the  depth  of  the  elastic  portion  in  compression,  b  and  d  are 
the  width  and  depth  of  the  section  respectively. 
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Equation  (3.1)  can  be  transformed  into  dimensionless 
form  by  dividing  by  the  yield  load  in  direct  compression,  Py 
which  is  equal  to  bdayc  to  give 

P/Py  =  1  -  0.5( 1-a/d) ( 1+aut/ayc )  (3.1) 

The  bending  moment  M,  on  the  cross-section  is  given  by 
M  =  JA  awdA 

=  0 . 5a y c ab ( d~a )  +  0 . 5ayc a '  (d/2-a-a ' /3 ) 
+0.5autb(d-a-a' )[d/2-(d~a-a' )/3] 

Simplifying  and  substituting  My  =  bd 2  cr  yc /6  ,  the  yield  moment 
in  pure  compression,  results  in 

M/My  =  (  1  -a/d )  (  1  +o  u  t  /or  y  c  )  (  0 . 5+a/d )  (3.2) 


Curvature 

Since  strains  are  small,  from  the  stress  diagram  of 
Fig.  2.2, 

tan  E0  —  E0  =  (Pyc+Pu  t  )  (3.3) 

(d-a ) 


The  limiting  curvature  for  elastic  conditions,  0y,  is 


0  y  =  My  =  2(7  y 
El  Ed 


(3.4) 
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From  Equations  (3.3)  and  (3.4) 

0  =  (  1+g|i  t/Pyr.  )  (3.5) 

0  y  (  1  -a/d) 


Substituting  Equations  (3.1)  and  (3.2)  in  (3.5)  gives  the 
following  M-P-0  relationship: 

0  =  4 ( 1 -P/P v ) 3  (3.6) 

0  y  [  3 (  1-P/Py  )-M/My ] 2 

Load  Curvature  Relationship 

Based  on  a  sine  curve,  the  equation  of  the  deflected 
shape  is 


y  =  6sin7rx/L 

where  y  is  the  deflection  at  any  given  distance  x  from  one 
end  of  the  member  of  length  L;  and  6  is  the  maximum 
deflection  at  midspan  as  shown  in  Figure  3.1.  The  curvature 
at  midspan  is  defined  as 

0  ( L/2  )  =  -y'  '  (L/2)  =  67T  2/L 2 

Using  the  expression  for  0y  given  in  equation  (3.4),  gives 

0m  =  tt2EI  6 
0 y  L2  My 

where  0m  is  the  maximum  curvature  occurring  at  column 
midspan.  Since  My =bd2 a y/6=Py d/6 ,  and  the  Euler  critical 
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load,  P  =7t2EI/L2  ,  then 


=  pr  65  (3.7) 

0  y  Py  d 

Equilibrium 

From  Equation  (3.6),  the  internal  moment  can  be  written 
as 

Mint  =  3(1-P/Py)  -  2 ( 1-P/Py )3/ 2  (3.8) 

My  (0/0 y ) ’ '  2 

The  external  moment  at  the  maximum  moment  section  (column 
midspan)  can  be  written  as 

Mext  =  M0  +P5  (3.9) 

My  My 

where  M0  is  the  moment  due  to  applied  loads.  Substituting 
Equation  (3.7)  into  Equation  (3.8),  and  equating  Mext=Mint 
yields 

M0  +  65  P  =  3 ( 1  —  P / P y )  2( 1-P/PV) 3/2  (3,10) 

d  Py  [ 6 ( 6/d ) (Pe/Py ) 3 1 7  2 


Equation  (3.10)  defines  the  failure  criterion. 


Expression  for  6/d 

Equating  Equations  (3.5)  and  (3.7),  for  0=0m  yields 
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(  1+gut  /o’  y  c  )  =  65  Pe  (3.11) 

2 ( 1 -a/d )  d  Py 

From  equation  (3.1) 

1-a/d  =  2 ( 1 -P/P v )  (3.12) 

(  1  +0  t  /<7y  ) 

Thus 

66  =  ( 1+q  t/av  ) 2PV  (3.13) 

d  4  (  1-P/Py  )  Pe" 

Substituting  Equation  (3.13)  into  Equation  (3.10)  and 
simplifying  yields  the  following  interaction  equation  at 
midspan 

M  =  3 (1-P/Py)  -  4 ( 1-P/PV ) 2  -  (1+q  t /a y  ) 2 P  (3.14) 

My  (1+q  t/oy  )  4(1-P/P')  Pe 

Again  from  equation  (3.1), 

q _ L  =  2 ( 1-P/Pv )  ~  1  (3.15) 

q  y  ( 1 -a/d) 

From  the  strain  distribution  of  Fig.  2.2 

a  =  6 u  c  (  1  -a )'  =  q y  c  (  l~q)  (  3  .  16) 

d  euc  +e  u<  a vc +a0 u t 

where  euc  is  the  ultimate  strain  in  compression  eut  is  the 
ultimate  tension  strain  and  a  is  the  ratio  of  yield  strain- 
to  ultimate  strain.  Substituting  Equation  (3.16)  into 
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Equation  (3.15)  and  setting  A=(1-P/Py)  and  A  =  aut/ayc 
a\2  +  2a( 1 -A ) X  +  (a-2A)  =  0 

from  which 

X  =  ( A— 1 )  ±  a2 ( A- 1 ) 2  +  a(2A-a)  (3.17) 

a 

Noting  that 

P  =  P/Pv 
P  P^TPy 

and  using  Py  =  bdae  E  and  P  =7r2EI/L2  and  I=bd3/12  then 

P  /Py  =  tt  2 /[  12(l/d)  2a6yc  ]  (3.18) 

Subtitution  of  equation  (3.18)  into  equation  (3.14)  yields 
the  final  interaction  equation  as 

M  =  3(1-P/Py)  -  4 ( 1 -P/Pv )  -  ( 1+gu t/a y  c ) x 1 2 P/P v ( L/d ) 2  ae vr 
My  (  1  +0  u  t/gy  c  )  4  (  1  -P/Py  )  7T2 

(3.19) 

Equations  (3.17)  and  (3.19)  can  be  used  to  generate 
interaction  diagrams  for  any  slenderness  ratio  and  material 
properties . 

3.3  Newmark's  Numerical  Integration  Procedure  (Method  2) 

This  method5,9  uses  a  more  precise  deflection  curve 
than  the  simple  sine  wave  assumed  in  Method  1.  It  is; 
however,  necessary  to  have  a  moment-thrust-curvature  (M-P-0) 
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relationship  for  the  type  of  material  and  cross-section 

being  analysed.  The  procedure  is  as  follows: 

1.  The  beam-column  is  divided  into  a  number  of  equal 
segments . 

2.  Values  of  L/d  and  P  are  selected  for  investigation. 

3.  A  low  value  of  M0  is  assumed. 

4.  A  trial  deflected  shape  is  assumed. 

5.  Using  P  from  Step  (2),  and  M0  from  Step  (3),  and  the 
deflections  from  Step  (4),  bending  moment  is  computed  at 
the  intermediate  points  in  the  span. 

6.  From  an  M-P-0  curve  of  the  material  and  cross-section, 
curvature  corresponding  to  the  total  moment  is  obtained. 

7.  A  new  deflected  shape  is  computed  by  integrating  the 
curvature  twice,  using  Newmark's  method  as  shown  in 
Figure  3.2. 

8.  Steps  (5),  (6)  and  (7)  are  repeated  with  the  new 

deflected  shape  until  convergence  to  a  fixed  shape  is 
obtained.  Fig.  3.1  shows  the  computation  procedure.  If 
M0  does  not  exceed  the  maximum  moment  the  member  can 
carry,  a  satisfactory  answer  is  obtained  after  3  or  4 
cycles . 

9.  Steps  (3)  through  (8)  are  repeated  with  a  new  value  of 
M0 

10.  End  slope  0o,  is  computed  for  each  M0 : 

@0  =  4W2/L  9 

where  w2  is  the  deflection  at  the  second  intermediate 
position  in  the  span  as  shown  in  Figure  3.2. 
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11.  The  M0-#0  curve  can  then  be  drawn;  or  the  highest  value 
of  M0  can  be  made  as  close  to,  but  yet  still  somewhat 
below  the  actual  value  of  (M0)max.  desired.  Fig.  3.3 
shows  a  typical  MO-0O  curve. 

The  above  procedure  was  applied  in  this  study  with  the 
following  peculiarities: 

1.  A  beam-column  with  four  point  lateral  loading  system  was 
analysed  as  shown  in  Figure  3.1. 

2.  The  beam-column  was  divided  into  eight  equal  segments. 

3.  The  deflected  shape  for  an  elastic  beam-column  with 
equal  end  moments  was  used  as  a  first  approximation  to 
the  true  deflected  shape. 

4.  The  axial  load  capacity  for  each  slenderness  ratio  was 
obtained  from  the  column  curve. 

5.  M-P-0  curves  obtained  from  material  properties  based  on 
small-scale  tests  were  employed,  as  shown  in  Figure  3.4. 

6.  The  criterion  for  maximum  moment  was  either  the 
exceeding  of  ultimate  moment  or  curvature,  or 
divergence  of  deflection  after  5  iterations. 

7.  The  final  deflected  shape  was  obtained  to  an  accuracy  of 
about  5%. 

3.4  Moment  Magnifier  Procedure  (Method  3) 

The  moment  magnifier  method5  is  an  approximate  equation 
for  the  design  of  beam-columns.  It  is  based  on  a  linear 
interaction  of  axial  load  and  moment  accounting  for 
additional  moment  due  to  axial  load  by  a  amplification 
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factor  1/(1-P/Pe).  It  is  adapted  to  predict  the  ultimate 
strength  of  timber  beam-column  in  the  following  form: 


£  +  CmM0  =  1  (3.20) 
Py  ( 1-P/Pe )My 

where  the  various  terms  are  as  defined  earlier,  and  Cm  is  a 
factor  dependent  on  the  distribution  of  the  applied  moment. 
For  the  type  of  loading  investigated  in  this  study  Cm  was 
taken  as  1.0,  as  it  was  close  to  a  uniformly  distributed 
loading . 

3.5  Computer  Codings  and  Interaction  Diagrams 

In  order  to  facilitate  the  tracing  of  the  interaction 
curves,  computer  programs  were  written  .for  each  of  the 
methods  described  above.  These  programs  are  placed  in 
Appendices  A1  through  A3,  corresponding  respectively  to 
Methods  1  through  3.  Various  values  of  material  properties 
and  slenderness  ratios  can  be  used  to  generate  as  many 
interaction  diagrams  as  desired. 

Interaction  curves  based  on  ultimate  compression  and 
tension  strains  of  0.0028mm/mm  and  0.0032mm/mm  respectively, 
a  value  of  a  of  0.864  and  slenderness  ratios  of  10,  20,  30, 
and  40  for  analysis  Method  2  are  shown  in  Figure  3.5. 
However,  the  material  may  possess  more  or  less  plasticity  in 
compression  and/or  a  different  tensile  strain  limit  at 
failure8.  Analysis  Method  2  is  therefore  examined  for  values 
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of  a  of  0.7,  0.8  and  0.9;  and  for  tensile  strains  of  30%, 
50%,  70%  and  100%  of  ultimate,  at  a  slenderness  ratio  of  20. 
Figure  3.6  indicate  reduction  in  moment  capacity  as  a 
decreases  from  0.7  to  0.9.  Figure  3.7  show  decrease  in 
moment  capacity  for  P/Py  less  than  0.3  as  the  tensile  strain 
at  failure  is  reduced  from  ultimate. 

3.6  Comparison  of  Analysis  Methods 

Figure  3.8  shows  interaction  curves  based  on  the  three 
analysis  methods  for  slenderness  ratios  of  10,  20  and  30.  It 
is  observed  from  the  curves  that  Method  1  gives  the  highest 
predictions.  Method  3  gives  the  least  moment  values.  Method 
2  seems  to  give  moment  values  intermediate  between  the  other 
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Figure  3.1  Computation  Procedure  for  Beam-Column  Deflection 
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Figure  3.2  Newmark's  Integration  Technique 
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Figure  3.3  Sample  Moment-Rotation  Curve 
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Figure  3.4  Moment-Curvature  Curves  for  Timber  Beam-Columns 
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Figure  3.5  Interaction  Curves  for  Timber  Beam-Columns 
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Figure  3.6  Interaction  Curves  for  Various  Depth  of  Yielding 

(Method  2) 
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Figure  3.7  Interaction  Curves  for  Various  Tensile 

(Method  2) 
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Figure  3.8  Interaction  Curves  Using  Various  Methods 


4.  EXPERIMENTAL  PROGRAM 


4.1  Test  Specimens 

4.1.1  Full  Scale  Specimens 

A  total  of  nine  full-scale,  factory  manufactured 
beam-column  specimens  were  tested.  The  test  specimens  were 
categorized  into  three  series  -  A,  B  and  C  corresponding  to 
three  cross-section  sizes:  175x152,  175x228  and  130x380.  The 
resulting  slenderness  ratios  (L/d)  were  13,  22,  and  33. 
Three  specimens  were  tested  in  each  series.  Table  4.1 
summarizes  the  properties  of  the  specimens. 

The  specimens  were  fabricated  in  a  plant  certified 
under  CSA  Standard  0177-M8113.  The  specimens  were  4990mm 
long.  The  laminations  were  38mm  thick  and  were  of  sufficient 
length  so  that  no  end  jointing  was  required.  Casco  1909 
cold-setting  casein  glue,  was  used  in  glueing  the 
laminations  together.  Shear  block  tests  were  performed  to 
confirm  the  adequacy  of  the  glue  bond. 

The  material  used  in  the  beam-columns  was  Douglas 
fir-Larch,  16c-E  grade  in  accordance  with  the  requirements 
of  CSA  Standard  0122-M8012.  The  average  modulus  of 
elasticity,  E,  for  all  laminations  equalled  or  exceeded  the 
minimum  of  12,400MPa  required  by  the  CSA  Standard.  Average 
moisture  content  for  all  laminations  ranged  from  7%  to  12%. 
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4.1.2  Small  Scale  Specimens 

Small  scale  specimens  for  establishing  tension  and 
compression  strengths  were  fabricated  from  material  cut  from 
three  pieces  selected  at  random  from  the  stock  used  for  the 
beam-column  specimens.  Figure  4.1  shows  the  cutting  pattern 
used  to  obtain  the  pieces  required  for  10  compression  and  10 
tension  specimens.  The  details  of  the  specimens,  based  on 
CSA  and  ASTM  standards  1 2 ' 1 4 ,  are  shown  in  Figures  4.2(a) 
and  (b).  The  moisture  content  for  all  the  small  scale 
specimens  ranged  between  6%  and  9%. 


4.2  Test  Set-Up 

The  beam-column  specimens  were  tested  in  a  horizontal 
test  frame  designed  to  provide  reactions  for  the  concentric 
axial  load  applied  through  a  hydraulic  jack,  rated  at 
4,450kN  maximum  load.  The  test  frame,  shown  in  Figures  4.3 
and  4.4,  consisted  of  rolled  wide  flange  sections  in  the 
longitudinal  direction  and  built-up  I-sections  as  cross 
members . 

Transverse  loads  were  applied  at  four  points, 
symmetrically  positioned  with  respect  to  the  midspan  of  the 
specimen.  The  lateral  loading  system  consisted  of  two  HSS 
sections  and  high  strength  steel  rods.  These  formed  a  yoke 
for  applying  pressure  through  load  cells  placed  at  positions 
PI  and  P2  shown  in  Plate  4.1.  Additional  HSS  sections  were 
used  as  distributing  beams  to  produce  a  four-point  loading 
system  on  the  specimen.  90kN  capacity  hydraulic  jacks, 
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reacting  against  the  laboratory  floor  slab  at  positions  PI 
and  P2  applied  loads  to  the  loading  yokes. 

Two  end  support  fixtures  held  the  ends  of  the  specimen 
in  position.  These  fixtures  transmitted  axial  load  to  the 
specimen  through  large  capacity  high  strength  steel  rotation 
rollers,  enclosed  by  machined  steel  plates.  Reaction  rollers 
placed  under  one  of  the  steel  plates  allowed  horizontal 
movement  of  the  end  of  the  specimen.  Because  of  the 
positioning  of  the  reaction  rollers,  temporary  adjustable 
supports  were  provided;  and  also  used  as  levelling  devices 
for  the  specimen. 

To  prevent  lateral-torsional  buckling  about  the  weak 
axis,  lateral  bracing  was  provided  close  to  positions  PI  and 
P2.  The  bracing  system  consisted  of  adjustable  HSS  sections 
acting  as  vertical  guides  to  the  specimen.  These  HSS 
sections  were  in  turn  braced  against  the  main  testing  frame. 
Figure  4.5  shows  a  diagram  of  the  loads  acting  on  the 
spec imen . 

4.3  Instrumentation 

Most  measurements  recorded  during  the  tests  were 
obtained  by  means  of  electronic  equipment.  Only  the  strain 
distribution  on  the  cross-section  was  measured  manually. 

An  electronic  load  cell  calibrated  to  a  maximum  load  of 
2,600kN  was  used  to  measure  the  axial  load.  160kN  capacity 
load  cells  were  used  for  measuring  lateral  loads.  The 
accuracy  of  the  load  measurement  is  considered  to  be  ±1%. 
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Deflections  were  measured  by  seven  electronic  linear, 
variable  differential  transformer  ( LVDT )  transducers.  A 
transducer  was  placed  at  each  of  the  four  lateral  load 
points,  one  was  placed  at  the  midspan,  and  the  remaining  two 
were  placed  as  close  as  possible  to  the  end  supports.  These 
end  transducers  were  used  to  monitor  the  effectiveness  of 
the  end  supports  and  provided  a  means  of  measuring  the  end 
rotations  by  measurement  of  deflections. 

Strains  at  midspan  were  obtained  by  means  of  a 
calibrated  125mm  Demec  gauge.  Demec  points  were  spaced  at 
25mm  for  the  A  and  B  series  and  at  40mm  for  the  C  series. 
Figure  4.6  shows  a  typical  arrangement.  Demec  gauge  readings 
were  recorded  manually.  All  the  data  from  the  load  cells  and 
transducers  were  fed  directly  into  the  NOVA  Computer  in  the 
laboratory . 


4.4  Test  Procedure 

The  specimen  was  positioned  in  the  end  support  fixture, 
and  filler  boards  were  placed  at  the  ends  to  ensure  full 
contact.  Due  attention  was  paid  to  proper  alignment  of  the 
specimen  to  ensure  concentric  axial  loading.  Wood  shims  were 
used  as  necessary  to  bring  the  specimen  to  the  required 
elevation . 

The  specimen  was  supported  temporarily  at  midspan  while 
the  vertical  loading  apparatus  was  positioned.  The  LVDT's 
were  positioned  and  Demec  points  were  installed. 
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To  ensure  that  the  loading  equipment  was  functioning 
properly,  pre-test  loads  of  approximately  10%  of  the  maximum 
axial  or  estimated  lateral  loads  were  applied.  The  pre-test 
loads  were  then  removed  and  any  necessary  adjustments  were 
made.  The  test  was  then  ready  to  begin. 

An  axial  load  of  about  20kN  was  applied  in  order  to 
hold  the  specimen  in  place  while  all  temporary  supports  and 
keeper  bars  were  removed.  A  set  of  readings  was  then  taken. 

The  axial  load  was  increased  in  increments  of 
approximately  one-fifth  of  the  total  axial  load,  with 
readings,  except  Demec  readings,  being  taken  at  each 
increment.  When  the  full  axial  load  was  reached,  it  was 
maintained  for  the  remainder  of  the  test.  A  complete  set  of 
readings  was  taken  at  this  point.  The  lateral  loads  were 
then  applied,  also  in  increments  of  one-fifth  their  expected 
maximum  value. 

At  each  increment,  the  loads  were  allowed  to  stabilize 
before  a  set  of  readings  was  taken.  However,  due  to  the 
nature  of  the  loading  equipment  (air  driven  motor  hydraulic 
pumps)  and  specimen  behaviour,  it  was  difficult  to  maintain 
the  load  at  a  precisely  constant  level. 

The  behaviour  of  the  specimen  was  monitored  by  plotting 
a  lateral  load  versus  midspan  deflection  curve  as  the  test 
progressed.  As  soon  as  the  specimen  began  to  show 
significant  non-linear  behaviour,  readings  were  taken  more 
frequently.  The  Demec  readings  were  discontinued  at  this 
stage.  All  specimens  were  tested  to  complete  failure. 
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The  test  data  were  transmitted  from  the  NOVA  computer 
to  the  AMDAHL  470  computer  for  further  processing. 


Table  4.1  Properties  of  Beam-Column  Specimens 
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(a)  Tension  Specimen 


610  mm 

< - «. 


(b)  Compression  Specimen 


Figure  4.2  Small-Scale  Specimens 
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Figure  4.3  Plan  View  of  Testing  Frame 
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Figure  4.4  Idealized  Test  Set-Up  and  Instrumentation 
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Figure  4.5  Beam-Column  Loading 
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Figure  4.6  Demec  Point  Spacing  on  Series  C  Specimens 
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Plate  4.1  Test  Set-Up 
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5.  TEST  RESULTS  AND  DISCUSSION 


5 . 1  Introduction 

Tables  5.1  and  5.2  summarize  the  results  of  compression 
and  tension  tests  on  the  small-scale  specimens.  Tables  5.3 
and  5.4  give  corresponding  summaries  for  the  full-scale 
tests . 

For  the  purpose  of  plotting  lateral  load  versus  midspan 
deflection  graphs,  the  average  of  the  loads  read  from 
load-cells  at  positions  PI  and  P2  was  used.  The  equipment 
load  amounted  to  a  total  of  1 . 1 kN  at  each  load  position.  The 
load-deflection  curves  for  all  specimens  are  shown  in 
Figures  5.1  to  5.3.  The  moment-end  rotation  curves  are  given 
in  Figures  5.4  to  5.6;  while  the  strain  distribution  on  the 
cross-sections  are  shown  in  Figures  5.7  to  5.13.  Figures 
5.15  to  5.22  show  the  deflected  shapes.  Figure  5.23  show 
test  results  together  with  analytical  predictions  for  33 
percent  ultimate  tension  strain,  using  method  2. 

During  the  tests,  the  laboratory  temperature  varied 
between  71  °F  and  74°F.  The  relative  humidity  varied  between 
26%  and  30%.  The  average  moisture  content  of  all  specimens 
was  approximately  7%. 

5.2  Small-Scale  Tests 

The  results  shown  in  Table  5.1  and  5.2  were  obtained 
from  the  10  compression  and  10  tension  tests  performed  on 
standard  specimens  as  described  in  Chapter  4. 
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In  compression,  the  average  yield  stress  was  48MPa. 
This  is  also  the  ultimate  strength,  using  an  elasto-plast ic 
stress  strain  curve.  The  average  ultimate  strain  was 
0 . 0028mm/mm ,  while  the  modulus  of  elasticity  averaged 


1 9 , 800MPa . 

The  value  of 

a  obtained  was 

0.864.  The  average 

values  in 

tension  were 

62MPa 

for 

ultimate 

stress , 

0 . 0032mm/mm 

for  ultimate 

strain 

and 

20 , 800MPa 

for  the 

modulus  of  elasticity. 

The  average  of  the  moduli  of  elasticity  in  tension  and 
compression,  20,300MPa,  was  used  in  all  computations. 
Coefficients  of  variation  for  all  calculated  averages  varied 
between  10%  and  22%. 

5.3  Full-Scale  Tests 

5.3.1  General  Behaviour 
Specimen  BCA 1 

Some  fine  cracks  were  observed  on  the  compression  face 
before  the  test.  These  cracks  did  not  seem  to  have  any 
significant  effect  on  the  behaviour  of  the  specimen. 

At  a  lateral  load  of  4kN,  a  cracking  sound  was  heard, 
at  a  load  of  7kN,  crack  openings  around  knots  and  knot-holes 
close  to  the  tension  face  initiated  failure.  At  a  load  of 
8kN,  splitting  in  the  tension  zone  was  becoming  common.  The 
beam-column  finally  failed  when  a  large  sloping  crack  formed 
through  the  bottom  two  laminations.  Plates  5.1(a)  and  (b) 
show  the  crack  patterns. 
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Specimen  BCA2 

When  the  lateral  load  due  to  the  equipment  was  applied, 
this  specimen  deflected  significantly  at  midspan.  At  a 
lateral  load  of  2kN,  a  wide  crack  developed  suddenly  from  a 
knot  in  the  second  lamination  from  the  bottom  face,  close  to 
load  position  PI.  This  crack  quickly  propagated  to  both 
sides  of  the  knot,  and  the  beam-column  became  very  sensitive 
to  the  lateral  load  adjustments.  An  attempt  to  increase  the 
load  caused  increasing  deflection,  resulting  in  final 
collapse  of  the  specimen.  Plates  5.2(a)  and  (b)  show  the 
observed  crack  patterns  at  failure. 

Specimen  BCA3 

A  number  of  knots  were  observed  throughout  this 
specimen.  As  for  BCA2,  significant  deflection  at  a  lateral 
load  of  1 kN  was  observed.  At  a  load  of  4kN,  cracks  started 
to  open  up  around  some  of  the  knots  close  to  midspan.  At  a 
load  of  6kN,  a  large  crack  opened  up  at  a  knot  in  the  bottom 
lamination  close  to  load  position  PI.  This  crack  penetrated 
three  bottom  laminations  along  a  sloping  grain,  to  cause  the 
final  failure.  Plates  5.3(a)  and  (b)  show  the  condition  of 
the  specimen  at  failure. 

Specimen  BCB1 

Initial  fine  cracks  were  observed  around  knots  at  a 
lateral  load  of  10kN,  close  to  load  position  PI.  Other  fine 
cracks  formed  at  a  load  of  15kN,  while  the  initial  cracks 
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opened  up.  Multiple  cracks  formed  at  midspan  and  propagated 
in  a  slightly  sloping  grain  to  cause  final  failure  at  a  load 
of  27kN.  Plates  5.4(a)  and  (b)  show  the  condition  of  the 
specimen  at  failure 

Specimen  BCB2 

The  behaviour  of  this  specimen  was  similar  to  that  of 
specimen  BCB1.  At  a  lateral  load  of  13kN,  a  knot  on  the 
compression  face  initiated  cracking.  The  specimen  failed  at 
a  load  of  14kN  when  large  inclined  cracks,  initiated  at 
small  knots  on  the  bottom  lamination,  opened  up 
significantly.  Plates  5.5(a)  and  (b)  show  the  final  failure 
condition . 

Specimen  BCB3 

This  specimen  behaved  similar  to  BCB1  and  BCB2.  Initial 
cracks  observed  prior  to  testing  on  one  side  of  the  specimen 
proved  of  no  significant  consequence.  At  a  load  of  6kN,  the 
beam-column  suddenly  cracked  around  two  knots  on  the 
compression  face.  As  a  load  of  7kN  was  reached,  compression 
failure  coupled  with  a  minor  edge  split  in  the  tension  zone 
caused  the  final  failure.  The  failure  cracks  are  shown  in 
Plates  5.6(a)  and  (b). 

Specimen  BCC1 

At  a  lateral  load  of  30kN,  minute  cracks  started  to 
open  up.  Compression  failure  occurred  close  to  load  position. 
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PI  at  a  lateral  load  of  40kN  and  near  load  position  P2  at  a 
load  of  45kN.  The  beam-column  failed  essentially  in 
compression,  with  minor  splitting  in  the  tension  zone  at  a 
load  of  54kN.  The  crack  patterns  are  shown  in  plates  5.7(a) 
and  (b) . 

Specimen  BCC2 

The  behaviour  of  this  specimen  was  somewhat  similar  to 
that  of  BCC1.  As  a  result  of  the  large  axial  load,  minute 
cracks  opened  up  before  application  of  any  lateral  load.  At 
a  lateral  load  of  17kN,  more  cracks  formed  in  the 
compression  zone.  These  cracks  widened  at  a  load  of  30kN.  At 
a  load  of  35kN,  the  beam-column  failed  by  crushing  in  the 
compression  zone,  confined  essentially  to  the  top  two 
laminations.  Plates  5.8(a)  and  (b)  show  the  crack  patterns. 

Specimen  BCC3 

This  specimen  carried  the  largest  axial  load,  which 
caused  small  cracks  to  open  up  at  zero  lateral  load.  At  a 
lateral  load  of  15kN,  the  cracks  widened  around  knots  close 
to  the  compression  face.  At  a  load  of  24kN,  crushing  in  the 
compression  zone,  coupled  with  some  tensile  cracks  at  a  knot 
caused  the  final  failure.  Plates  5.9(a)  and  (b)  show  the 
conditions  at  failure. 
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5.3.2  Load  Deflection  Curves 
Series  A  Specimens 

The  load-deflection  curves  obtained  for  the  series  A 
specimens  are  shown  in  Fig.  5.1.  As  was  expected,  the 
lateral  load  at  failure  decreased  as  the  concentric  axial 
load  was  increased.  A  general  non-linear  but  relatively 
smooth  plot  for  specimen  BCA 1  was  perhaps  due  to  the  light 
axial  load.  The  behaviour  of  specimen  BCA3  was  irregular  due 
to  the  significant  number  of  knots  present,  and  perhaps  the 
higher  level  of  axial  load.  Specimen  BCA2  with  the  largest 
axial  load  in  this  series  (67  percent  of  the  Euler  load), 
failed  in  an  instability  mode  as  indicated  by  the  load 
deflection  curve. 


Series  B  Specimens 

The  behaviour  of  Series  B  specimens  was  similar  to  that 
observed  for  the  series  A  specimens.  The  stockier  nature  of 
these  specimens  is  however  reflected  in  the  relatively 
smoother  load-deflection  curves  as  shown  in  Fig.  5.2.  For 
specimen  BCB3,  significant  deflection  occurred  when  the 
lateral  load  due  to  the  equipment  was  applied.  Also  at  a 
lateral  load  of  approximately  4kN,  an  apparent  defect  caused 
a  significant  increase  in  deflection. 


Series  C  Specimens 

Series  C  specimens  which  were  the  stockiest  specimens 
produced  the  smoothest  load-deflection  curves  as  shown  in 
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Fig.  5.3.  However,  the  effects  of  knots  and  knotholes  were 
still  evident.  Specimen  BCC3  showed  some  upward  deflection 
after  the  introduction  of  full  axial  load.  This,  however, 
did  not  significantly  affect  its  final  failure. 

5.3.3  Moment-End  Rotation  Curves 

Figures  5.4  to  5.6  show  the  moment  versus  end  rotation 
curves  for  all  specimens.  As  expected  the  shapes  of  these 
curves  are  similar  to  those  of  the  load-deflection  curves. 

5.3.4  Strain  Distribution 

The  measured  strains  across  the  cross-section  for 
specimen  BCA1,  series  B  and  C  specimens  are  shown  in  Figures 
5.7  to  5.13.  It  was  difficult  to  measure  strains  on  the 
cross-section  of  specimens  BCA2  and  BCA3  because  of  unstable 
behaviour.  The  effect  of  knots  and  knotholes  on  the  linear 
distribution  of  strain  could  be  observed  in  some  of  the 
figures . 

5.3.5  Deflected  Shapes 

Figures  5.15  to  5.22  show  the  deflected  shape  for  all 
specimens.  It  is  observed  that  these  shapes  are  similar  to  a 
portion  of  sine  wave  in  most  cases,  especially  at  loads 
close  to  failure.  Effects  of  knots  or  knotholes  sometimes 
affected  the  shapes  as  evident  in  some  of  the  figures. 
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5.4  Typical  Failure  Modes 

Tensile  splitting  produced  by  sloping  grain  was  the 
common  failure  mode  of  the  series  A  specimens.  This  was 
perhaps  due  to  the  high  slenderness  ratio  and  presence  of 
knots  and  knotholes.  Strain  measurements  on  the 
cross-section  of  specimen  BCA1  indicated  compression  failure 
preceding  tensile  splitting.  Because  of  unstable  behaviour 
of  specimens  BCA2  and  BCA3 ,  this  failure  mode  could  not  be 
conf i rmed . 

Specimens  BCB1  and  BCB2  were  the  best  examples  of  the 
assumed  failure  criterion  employed  in  the  analyses.  Crushing 
close  to  the  compression  face  followed  by  tensile  splitting 
was  observed  for  both  specimens. 

Specimens  BCB3  and  all  series  C  specimens  failed 
essentially  in  compression,  usually  around  knots  and 
knotholes.  However,  it  was  interesting  to  observe  slight 
splitting  in  the  tension  zone,  except  for  specimen  BCC2 
which  showed  no  sign  of  any  tensile  splitting  (See  Plate 
5.8).  The  cross-section  of  specimen  BCC3  was  in  compression 
up  to  70%  of  its  failure  load,  as  indicated  in  Fig.  5.13. 

5.5  Ultimate  Strength  of  Beam-Columns 

Table  5.3  shows  the  maximum  strengths  of  the 
beam-column  specimens  based  on  properties  listed  in  Table 
4.1.  The  mean  value  of  the  axial  force  (Col.  2),  which  was 
maintained  practically  constant  during  a  given  test,  was 
used  in  estimating  the  bending  moment  due  to  axial  load. 
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(Col.  7).  The  maximum  value  of  the  total  bending  moment 
(Col.  8)  is  then  the  sum  of  the  moment  due  to  transverse 
load  (Col.  5)  and  the  moment  from  Col.  7. 

5.6  Comparison  of  Test  Results  With  Analytical  Predictions 

Table  5.4  shows  the  measured  and  predicted  strength  for 
all  beam-columns,  using  all  analyses  Methods  (Cols.  6  and 
7).  The  predictions  based  on  Method  1  (Col.  3)  are  observed 
to  be  generally  within  10%  of  those  based  on  Method  2  (Col. 
4).  This  is  not  surprising  since  for  most  specimens, 
measured  deflected  shapes  were  similar  to  the  sine  waveform 
assumed  in  analysis  Method  1.  It  is  also  observed  that 
analysis  Method  3  (Col.  3)  gives  predicted  moments  close  but 
generally  conservative  compared  to  those  from  Method  2.  It 
thus  has  potential  for  use  as  a  method  of  design  because  of 
its  simplicity.  The  result  for  specimen  BCA2  is  questioned 
because  of  its  instability  failure.  Also,  Specimen  BCB3 
failed  prematurely  in  compression,  limiting  its  capacity 
compared  to  the  other  two  Series  B  specimens.  This  result  is 
also  questioned.  The  difference  between  predicted  and  test 
moments  may  be  attributed  to: 

1.  The  influence  on  ultimate  strength  of  natural  defects 
such  as  knots,  knotholes  and  sloping  grain.  This  is 
evident  from  the  strain  distribution  on  the 
cross-section  of  the  various  specimens.  Only  Series  B 
specimens  attained  an  average  of  about  40%  ultimate 
tension  strain  at  failure. 
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2.  The  usual  scatter  in  results  from  timber  tests,  which 
does  not  allow  for  much  confidence  in  the  use  of  average 
values.  It  is  observed  that  the  average  ultimate  strains 
used  in  the  prediction  curves  have  a  coefficient  of 
variation  of  more  than  20%.  Perhaps  matched  small-scale 
specimens,  from  an  increased  number  of  full-scale  tests, 
may  give  closer  agreement  between  analysis  and  test 
results . 

3.  The  relatively  small  number  of  tests  reported  in  this 
investigation . 

To  account  for  the  above  effects  the  analyses  may  be 
modified  by  applying  an  undercapacity  factor.  To  illustrate 
this  approach  an  undercapacity  factor  of  0.7  has  been 
applied  to  analysis  Method  2  (Col.  6  of  Table  5.4).  The 
resulting  modified  interaction  diagrams  are  shown  in  Figure 
5.23  together  with  test  results.  With  the  exception  of  the 
questioned  test  results,  a  good  correlation  is  observed. 
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Table  5.1  Compression  Test  Results 
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Table  5.2  Tension  Test  Results 
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Table  5.3  Beam-Column  Test  Results 
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Table  5.4  Measured  and  Predicted  Values 
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Figure  5.1  Load-Deflection  Curves  for  Series  A  Specimens 
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Figure  5.2  Load-Deflection  Curves  for  Series  B  Specimens 
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Figure  5.3  Load-Deflection  Curves  for  Series  C  Specimens 
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Figure  5.4  Moment-Rotation  Curves  for  Series  A  Specimens 
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Figure  5.5  Moment-Rotation  Curves  for  Series  B  Specimens 
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Figure  5.6  Moment-Rotation  Curves  for  Series  C  Specimens 
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Figure  5.7  Strain  Distribution  on  Specimen  BCA 1 
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Figure  5.8  Strain  Distribution  on  Specimen  BCB1 
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Figure  5.9  Strain  Distribution  on  Specimen  BCB2 


DEPTH(mfn) 


66 


STRRIN(mm/mm)  xiO 


4000 


Figure  5.10  Strain  Distribution  on  Specimen  BCB3 
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Figure  5.11  Strain  Distribution  on  Specimen  BCC 1 
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Figure  5.12  Strain  Distribution  on  Specimen  BCC2 
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Figure  5.13  Strain  Distribution  on  Specimen  BCC3 
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Figure  5.14  Deflected  Shape  of  Specimen  BCA 1 
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Figure  5.17  Deflected  Shape  of  Specimen  BCB1 
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Figure  5.20  Deflected  Shape  of  Specimen  BCC 1 
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Figure  5.23  Interaction  Curves  Based  on  Undercapacity  Factor 

of  0.7 
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Plate  5.1  Failure  of  Specimen  BCA 1 
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Plate  5.2  Failure  of  Specimen  BCA2 
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Plate  5.3  Failure  of  Specimen  BCA3 


. 


: 


(a ) 


(b) 


Plate  5.4  Failure  of  Specimen  BCB1 
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Plate  5.5  Failure  of  Specimen  BCB2 
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(b) 


Plate  5.6  Failure  of  Specimen  BCB3 
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Plate  5.7  Failure  of  Specimen  BCC 1 
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Plate  5.8  Failure  of  Specimen  BCC2 
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Plate  5.9  Failure  of  Specimen  BCC3 
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6.  CONCLUSIONS  AND  RECOMMENDATIONS 


6.1  Observations  and  Conclusions 

The  following  are  the  results  of  the  theoretical  and 

experimental  investigation  of  timber  beam-columns  reported 

herein : 

1.  The  theoretical  interaction  curves  relating  the  bending 
momentfM,  and  the  compressive  axial  force,  P,  depend  on 
the  amount  of  plasticity  in  the  compression  zone  and  the 
tension  strain  at  failure. 

2.  The  theoretical  prediction  based  on  sine  wave  deflected 

shape  of  beam-column  gives  the  highest  values  of 

moments.  The  moment  magnifier  approach  gives  the 

smallest  moment  values.  The  Newmark's  numerical 
integration  method  gives  moment  values  intermediate 
between  the  other  two  methods. 

3.  The  value  of  the  tensile  strain  at  failure  is  dependent 
on  the  type  of  material  and  the  magnitude  of  the  applied 
axial  load. 

4.  The  interaction  curves  obtained  from  the  Newmark's 
numerical  analysis  pocedure  and  modified  by  an 
undercapacity  factor  of  0.7  are  in  good  correlation  with 
the  test  results. 

5.  The  moment  magnifier  approach  has  potential  as  a  design 
method  if  used  with  an  appropriate  undercapacity  factor. 

6.  Insufficient  specimens  were  tested  to  provide 
statistically  significant  data.  The  present  program  must 
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be  viewed  as  a  pilot  study. 


6.2  Recommendations  for  Further  Study 

The  following  recommendations  should  be  considered  in 
further  investigations  into  the  ultimate  strength  behaviour 
of  timber  beam-columns: 

1 .  A  number  of  tests  should  be  conducted  at  each  axial  load 
level  for  each  slenderness  ratio. 

2.  A  larger  number  of  small-scale  specimens,  more  closely 
related  to  each  beam-column  specimen  should  be  tested  in 
order  to  better  define  the  material  properties  of  each 
beam-column  specimen. 

3.  During  fabrication,  care  should  be  taken  to  ensure  that 
obvious  defects  are  not  located  at  critical  sections  and 
in  critical  laminations. 

4.  The  use  of  electrical  resistance  strain  gauges  to 
measure  compression  and  tension  strains  should  be 
considered. 
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4 
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7 
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9 

1  O 

1  1 

1  2 

1  3 

1  4 

1  5 

1  6 

1  7 

1  8 

1  9 

20 

2  1 

r  i 

23 

2  4 

25 

2  6 

27 

2  £ 

29 

30 

3  1 

32 

33 

34 

35 

36 

37 

3  £ 

3  S 

40 

4  1 

«  2 

43 

4  4 

4  5 

4  6 

4  7 

4  £ 

4  9 

5  C 

5  1 

52 

53 

54 

5  5 

56 

5  7 

5  £ 

5  9 

6  0 

6  1 

E  2 

E  2 

6  4 

6  5 

6  E 

6  7 

6  6 

6  9 

70 

7  1 

72 

73 

7  4 

75 

76 

77 

78 

79 

80 

8  1 

8  2 

8  3 

6  4 

£  5 

8  6 

8  7 

88 

8  9 

9  0 

9  1 

S  2 
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C  THIS  PROGRAMME  COMPUTES  INTERACTION  DIAGRAM 

C  FOR  WOOD  BEAM-COLUMNS  GIVEN  VARIOUS  MATERIAL 

C  PROPERTIES  USING  A  SINGLE  EQUATION 

C  IT  ALSO  PLOTS  THE  OUTPUT 

C**»sssss**sss«s*«,,,,*s. «*«**, ssx, 

DIMENSION  PO VP Y ( 40 ) , EMO VMY < 40 ) 

C 

C..  PRINT  TITLE  AND  READ  INPUTS 
WRITE (6 , 2000 > 

J=  1 

1  R  E  A  D  <  5  ,  1000.END=999  I  ELOVD , ALPHA ,  EUC.EUT  ,  FR AC 
M=  C 

EYC=ALPHA*EUC 
ET  sFRAC*EUT 
SRL IMsET/EYC 
C 

C  PRINT  PARAMETERS  AND  TABLE  HEADS 

WR  I  T  E ( 6 , 21 OOIELOVD, ALPHA,  EUC.ET.SRLIM 
WR  I TE  <  6  , 2200  > 

C 

C.  COMPUTE  THE  DIFFERENT  LOADS 
IF(ELOVD.LE.O.  )  GO  TO  1 
P  1  =  3  .  14  15  9  3 
P I S0=  P I  *  *  2 

N  D  s  1 

C 

C  COMPUTE  P/P Y , S  I  CT/S I C  .  UC  , MO/MY 
10  DO  100  1*5.105,5 

POVPY  ( NO  )  =  F  L  0  A  T (  1-5  1/100 

IF  *  POVPY  (NO  >  EO  .  1  .  )  POVPY  ( ND  1  =0 . 999 

A=  (  1  .  -POVPY  (  ND  )  ) 

BE  =  2  « A* ALPHA  -  2  .* ALPHA 
DlSCsBE**2*<  »ALPHA» < 2 . » A- ALPHA ) 

I  F  (  D I  SC  .  LT  .  O  .  1  GO  TO  100 
DErSORT  (DISC) 

SR*  (  B  E  +  D  E  )  /  (  2  .  ‘ALPHA ) 

IF(SR . GT ,SRL1M)SR=SRLIM 
C  .  .  CALCULATE  MO/MY 

EMOVMY  (  ND  >  *3  «A-4  *A“2/(1  *  S  R  >  -  (  1  ♦  S  R  '  *  *  2  /  <  4  «£) 

**POVPY<  ND  >  *  12  .  *EUC*ALPHA»  <EL0VD“2  >/P!SC 

C  .  OUTPUTTHE  RESULTS 

I F  (  EMO VMY  (  ND  )  .  L T  .  O  .  > GC  TO  15C 
WR  I  T  E  (  6  ,  2300  '  POVPY  (  ND  •  , SR  ,  EMO VMY  l  NO  I 
ND*  ND*  1 
100  CONTINUE 
CO  TO  2  O  C 

120  IFIPOVPYINDI  LT.O  (GO  TO  998 
150  POVPY ( ND » =POVPY( ND > -O  OO 1 
A*  1  - POVP Y  <  ND  » 

BEs 2  * A*ALPHA - 2  ‘ALPHA 
DlSC*BE*»2  +  4  .  sALPHA*  l  2  *  A •  A  L  P  H  A  ) 

IFID1SC  .LT.O.  )  GP  TC  99S 
DE=SORT (DISC) 

SRs  (BE  +  DE  )  /  ( 2  .  ‘ALPHA) 

IF (SR . GT .SRLIM)SR=SRL1M 
C  .CALCULATE  MO/MY 

EMOVMY  <  ND  )*3  *  A  -  4  *  A *  *  2 / (  1  .  ♦ S R  >  *  (  1  +SR‘**2/<4  *A> 

«*POVRY ( NO  )*  12  *EUC‘ALPHA*  (  E L 0 VO* *  2  )  / P  I S C 
IF  (  POVPY  <  NO  )  .  LT  . 0 .  > GO  TO  200 
C 

C  .  OUTPUT  THE  RESULTS 

160  WP.  ITE  (  6 , 2300  IPOVPY  (ND  '  ,  SR  ,  EMOVMY  (  ND  ' 


2  OC  CALL  PLOTIT  ( EMOVMY  ,  POVPY  ,  ND  ,  J ,  1  ,  1  .3,0  0,0.25. 
*5  .0.  .0.15,8  6  I 

J=  J+  1 

C  GC  TO  READ  FRESH  INPUTS 
GO  TO  1 

998  WR  7  T  E  (  6  ,  2 400  » PO VP Y I  ND  >  ,  EMOVMY  (  ND  ' 

9  9  9  Js  O 

CALL  PLOTIT (EMOVMY , POVPY , ND , J, 1 , 1 ,3,0.0,025, 
*5  C  ,0.15,8  .6) 

STOP 


C 

C  .FORMAT  STATEMENTS 
1 OOO  FORMAT  (  5  F  1 0  .  S  . 

2000  F0RMAT(///,23X, 'INTERACTION  CALCULATIONS  FOR  WOOD'// 
*23X , 'BEAM-COLUMNS-  RECTANGULAR  SECTIONS'// 

*  3  4  X  ,  'WI NTER  1  A  '  / / 

*  2  3  X  ,  '****»»**»«****»***»******““““*'  ) 

2100  FORMAT  (/// 10X  ,'  L/D*  '  ,  F  6  .  5  ,  5  X  ,  'ALPHA:  '  ,  F 8  .  5  ,  / /  1 O X  ,  '  E U C  = 

*  F  £  5  ,  5  X  ,  '  E  T  =  '  ,  F  8  5,5X,  'SRLIM:  ',F£  5//) 

2200  FORMAT  (  10X.  '  P/P Y  '  ,  10X.  '  S  I  G  .  T /S 1 C  . UC '  ,9X,  'MO/MY'//  I 
2300  FORMAT  (  1  OX  .  F5  .  3  ,  1 2X  .  F6  4  ,  1  1 X  ,  F7  .  4 / /  ' 

2400  F0RMAT(//1OX,  'NEGATIVE  VALUE  OF  P/P Y  *',F10.6, 
s'AND  M/MY  s  '  ,  F  1  O  6  /  ‘ 

STOP 

END 


file 


- 


. 


■ 

■ 
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Appendix  A2 


Method  2 


1 

2 

3 

a 

5 

6 

7 

6 

9 

1  0 

1  1 

1  2 

1  3 

1  4 

1  5 

1  E 

1  7 

1  6 

1  9 

2  C 

2  1 

22 

23 

24 

25 

26 

27 

26 

29 

30 

3  1 

32 

33 

34 

35 

36 

37 

36 

39 

40 

4  1 

4  2 

4  3 

4  4 

4  5 

4$ 

4  7 

46 

4  S 

50 

5  1 

52 

53 

54 

55 

56 

57 

5  6 

5  S 

6  0 

6  1 

62 

6  3 

64 

6  5 

6  6 

6  7 

6  £ 

6  9 

70 

7  1 

72 

73 

7  4 

75 

76 

77 

78 

7  9 

80 

8  1 

82 

63 

84 

6  5 

8  6 

87 

86 

89 

90 

9  1 

9  2 

9  3 

9  4 

9  5 

9  6 

9  7 

9  6 

9  9 

1  00 

1  O  1 

102 

103 

104 

1  05 

1  06 

1  07 

1  06 

1  09 

1  1  0 

1  1  1 

1  1  2 

1  1  J 
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C**i 

c  . 

c 

c 

c 

c 

c 

c 

c 

c 

c  *  *  * 
c 


THIS  PROGRAMME  CALCULATES  THE  ULTIMATE  STRENGTH  OF  TIMBER 
BEAM-COLUMNS  BY  NEWMARK'S  PROCEDURE  IT  IS  THE  LATEST  VERSION 
OF  W1NTER4.  MODIFIED  FOR  A  FOUR  POINT  LATERAL  LOADING  SYSTEM 
IT  CALCULATES  THE  FINAL  DEFLECTION  TC  5  PERCENT  ACCURACY. 


INTEGER  C , 0 , V 

COMMON  W(  10)  ,  WO  <  10)  ,  EMOM Y  <  10)  ,EMOMYT(  10)  ,PDEL  (  10)  .PHOPHY  I  10)  , 
*SLOPE(  lO)  ,Ri  10)  , EY  ,  EL  OVD  ,  POVPY  ,  A  ,  PHUPHY  ,  EOOOY  ,  EMOPY  (  10)  , 

*  EM  0  V  P  (  10)  ,  V  ,  EMYOPY  ,  NN  ,  EMUT  ,  EMOMYM (  40)  ,THETAO(40)  .POVPYM.JJ 
WR  I  TE  <  6  .  2000  ) 

READ (5,950 MRU),  1  =  1,9) 

5  READ(5,9O0,END=999)EL0VD,ALPHA,EUC,EUT,FRAC 
EYs ALPHA* EUC 
ET  =  FR AC  »  E  UT 
SRL IM=ET /EUC 

WRITE(6,21OO)EL0VD, ALPHA, EY. EUC. ET.SRLIM 

..CALCULATE  PURE  AXIAL  CAPACITY  FROM  COLN.  CURVE 
CALL  STAR 
CALL  AXCAP 
CALL  STAR 

CALL  PLOT  IT  ( THETAO ,  EMOMYM .  ND  .  NP  ,  1  ,  1  ,  3 , 0  .  .  ,5  ,0..  ,6  ,6) 

10  READ < 5  ,  1 OOO  ,  END  =  99  9  )  POVPY 
..CHECK  FOR  ANOTHER  L/D  FOR  CONSIDERATION  -  NEC  P/P Y 

I  F  (  POVPY  .  LT  .  O  )  GO  TO  5 

IS  NP  =  NP+  1 
0  P  C  T  =  0 
A  *  1  .  -  POVPY 

BE=2  *A*ALPHA-2  ‘ALPHA 
D I SC  =  BE *  *  2  +  4  ‘ALPHA*  (2  ‘A-ALPHA  ) 

IFIDISC.LT.O. ) WR I TE ( 6 , 2600 ' 

DE=S0RT(DISCI 

SR = ( BE +DE ) / ( 2  ‘ALPHA! 

SRL IM=  ET /EY 

I F  ( S  R  .CT  SRLIM)SR=SRLIM 
A  0  V  D  =  (  EUC*  <  1  .  -ALPHA)  )  /  <  EUC  +  ET  ) 

PHUPHYs  !  1  .  ♦ S  R  )  /  ( 2  *  <  1  .  -AOVD)  ) 

EMUTs  3  .  *  A  -  (2  *  ( A  >  **  1  .  5  I/SORT  (  PHUPHY  ! 

EMYOPYs  1  .  /  <  6  ‘ELOVD  ) 

.INITIALISE  VARIOUS  COUNTERS 
L  =  1 
M=  O 
C  =  O 
V  =  o 
NPsO 

II  =  1 
J  J=  1 
N  D  =  0 

CALL  STAR 
WR  I  T  E (  6 ,2150; 

CALCULATE  ASSUMED  MOMENT 
20  CALL  MOMENT 

.CHECK  IF  IN  FINER  M/MY  SECTION 
IF(C.CT.O)GO  TC  550 
30  CALL  ELADEF 

CALCULATE  THE  PDELTA  EFFECT  AND  ADD  THE  EFFECT 

50  CO  70  1=1,9 

PDEL ( I ) sPOVPY * WO ( I ) /EMYOPY 
EMOMYT  (  I  ) =  EM  OM  Y  (  I  >  ■*•  P  D  E  L  (  I  ) 

70  CONTINUE 

.CALCULATE  PHI  AND  DEFLECTION  FROM  M-P-0  RELATIONS 
CALL  PHIDEF 

CORRECT  THE  DEFLECTIONS 
CALL  CORDEF 
Lb  l*  1 

CHECK  IF  CONVERGENCE 

T  0  L  ER  =  0 
DO  100  J=  2  .  8 
I F ( WO l J ) . EO . O ) GO  TO  1 OO 
D I F  =  ABS ( W  <  J  ►  -WO l J  )  ) 

OPCT=D I F/WO ( J ) 

100  TOLER  =  TOLER+ABS(OPCT  J 

I  F ( TO l ER  .  L E  . O  . 2 1 O  OR  .  L  .  CT  .  5  ) GO  TO  300 


.CHECK  IF  MAXIMUM  CURVATURE  EXCEEDED 
1  F  <  PHOPHY  <  NN  ).  GT  .  PHUPHY  )  GO  TO  500 


OTHERWISE  USE  NEW  DEFLECTION  ESTIMATE 


DO  200  K= 1 , 9 
200  WO ( K ) =W( K ) 


. 

. 

- 


I 

1 


- 


1  1  4 

1  1  5 

1  1  6 

1  1  7 

1  1  8 

1  1  9 

1  20 

1  2  1 

1  22 

1  23 

1  24 

1  25 

1  26 

1  27 

1  26 

1  29 

1  30 

1  3  1 

132 

1  33 

1  34 

1  35 

1  36 

137 

1  36 

1  3  9 

1  40 

1  4  1 

1  42 

1  43 

1  4  4 

1  4  5 

1  46 

1  4  7 

1  4  8 

1  4  9 

1  5  C 

1  5  1 

1  52 

1  53 

1  5  < 

1  55 

1  56 

1  57 

1  56 

1  59 

1  6  0 

1  E  1 

1  6  2 

1  63 

1  6  4 

1  6  5 

1  6  6 

1  6  7 

1  6  £ 

1  6  9 

1  70 

1  7  1 

172 

1  73 

1  74 

1  75 

1  76 

1  77 

1  7  £ 

1  75 

1  60 

1  6  1 

1  82 

1  8  3 

1  84 

1  85 

1  66 

1  87 

1  86 

1  85 

1  90 

1  9  1 

1  92 

1  9  3 

1  9  4 

1  95 

1  9  6 

1  9  7 

1  9  £ 

1  9  9 

200 

20  1 

202 

203 

204 

205 

206 

207 

206 

209 

2  1  0 

2  1  1 

2  1  2 

2  1  3 

2  1  4 

2  1  5 

2  1  6 

2  1  7 

2  1  8 

2  1  9 

220 

22  1 

2  2  2 

222 

224 

225 

2  2  S 
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OPCTr O 
CO  TO  50 

300  IF  < TO l ER . I E . O . 2 1 O > GO  TO  400 
CALL  STAR 

WRITE<6 , 2  200  )  EMOM Y  (  5  ) 

CALL  STAR 
C 

C  CHECK  IF  MAX.  M/MY ( AT  POINT  5>  IS  LESS  THAN  0  10 

I F ( EMOM Y ( 5 ) . L T . O  1)G0  TO  10 
1FIC.GT.0)  GO  TO  10 
GO  TO  500 

C  CALCULATE  THETAO  FOR  CONVERGENCE  ACCORDING  TO  CHEN  AND  ATSUTA 
400  THETAOI 1 ] 1*4  *W(3 I 

N  I  T«L -  1 

EMOM  YM (  1  I  ) sEMOMY  <  5  ) 

WR  I  TE (  6 , 2300  ) PO VP  Y  ,NIT  ,  EMOM YM (  I  I  )  ,  THETAO  <  I  1  )  ,  TOLER 
C 

C.. CHECK  FOR  NEGATIVE  THETAO 
C 

IF  ( THETAO (  I  I  )  .  L T  .  O  .  )  CO  TO  10 
C 

4  50  N  D  =  N  D  + 1 
11*114-1 

IFIEM0MYI5  i  .  CE  .  EMUT)WRITE<6 ,2500)  EM UT 
I F ( EMOMY < 5 ) . GT . EMUT ) CO  TO  10 
C 

C  REPEAT  CALCULATION  FOR  OTHER  M/MY  VALUES 
C 

IFIC.GT.OIGO  TO  550 
C  REINITIALISE  COUNTER 
L  *  1 

0PCT=O 
CO  TO  20 

c 

c  .USE  FINER  FRACTIONS  OF  M/MY  BEFORE  DIVERGENCE 
C  REINITIALISE  COUNTER  AGAIN 
C  CHECK  IF  MAX  M/MY  IS  LESS  THAN  0.1 

c 

500  1FIEM0MYI5 I ,LT .0. 1 IGO  TO  10 
IFIC.GT.OIGO  TO  550 
L  *  1 
C 

C  TAKE  BACK  MOMENT  CALCULATION  ONE  STEP 

c 

V  =  V  -  4 
C  =  C  ♦  1 
CO  TO  20 

550  DO  600  0=1, S 

EMOMY  (O  )  =  EMOMY  (0)4-0 .01 
I F ( EMOM Y ( O > . L T  O.)C0  TO  10 
600  CONTINUE 
C 

C  .CONITNUE  FOR  ANOTHER  SET  OF  ITERATIONS  IN  FINER  FRACTIONS 


C  =  C  4-  1 

IFlC.GE  1OIC0  TO 
0  P  C  T  =  C 
GO  TO  30 
C 

c.  FORMAT  STATEMENTS 
C 

900  FORMAT (5FI0.5I 
9  SC  FORMAT  (  1  OF  £  6  » 

1  000  FORMAT  ( F  1 C  .  5  » 
2000  FORMAT < 25X , 'WOOD 
* ' NEWMARKS  METHOD 


1  0 


ULTIMATE  STRENGTH 
WINTER4C' ,/// » 


BY', 


2100  FORMAT  I//10X,  'L/D=  '  ,F10  3,5X,'ALPHA=  '  ,  F 1 0 . 5  ,  5 X  ,  '  E Y = 

*  F 1 O  5.//10X,  '  E  U  C  =  '.F10.5.5X,  'ETs  '  , F 1 0  .  5 , 5X  ,  ' SR L 1 M=  '  ,  F  1  C  5  ) 

2150  FORMAT  (  / / 1  OX  ,  'P/PY'  .  lOX,  'NIT'  ,  10X,  'M/MY'  .lOX,  'THETAO'  , 

* 10X . 'TOLERANCE ' ) 

2200  FORMAT!  12X,  'DIVERGENCE  OF  DEFLS  AFTER  5  ITERATIONS', 

*  ' AT  MAX  .  M/MY  =  '  ,  F7  4  i 

2300  FORMAT  (  / 1  OX  ,  F4  3  ,  1  OX  ,  I  3  ,  9 X  ,  F5  .  3  ,  8X  ,  F 1 O  8  ,  9 X ,  F 6  4) 

2500  FORMAT ( /'*i******i*MO  IS  GREATER  THAN  ULTIMATE, MUL  =  ' ,F10  7) 

2600  FORMAT (/'********** *NEGAT I VE  DISC  IN  SIGT/SIGUC') 

999  STOP 
END 
C 

C  .  .  SUBROUTINE  TO  CALCULATE  PURE  AXIAL  CAPACITY  FROM  COLN  CURVE 

Czxzzzzzzzzzzxzzzzzxzxxxzzxzxzxxzzzzzxzzxzzzzzxzzzzzzzzzzzzzzzzz 

SUBROUTINE  AXCA» 

I  NTEGER  C  ,  0 .  V 

COMMON  W<  lO)  ,W0<  10)  ,  EMOMY  (  10)  ,  EMOMY  T  <  10)  ,  PDEL  I  10)  .PHOPHYI  10)  , 

*  S  L  0  P  E  I  10)  ,  R  (  10)  ,  EY  ,  ELOVD,  PO VP  Y , A  ,  PHUPHY  ,  EOOOY  ,  EMOPY  (  10)  . 

»  EM 0  V  P  (  1O/,V,EMY0PY,NN,EMUT,EM0MYM(4O)  ,  THETAOUOl  ,  POVPYM,  JJ 
C 

C  .  ENTER  L/D  VALUE  SEPARATING  SHORT  FROM  LONG  COLUMN 
C..  NOTE  THE  ASUMPTION  OF  SIGUC=SIGY  IN  CALCULATION 
C 

SUOVSYsl . OO 
Kt  2  1 

c 

C  .CHECK  SLENDERNESS  RANGE 
C 

IF(ELOVD.LE  K  )  G  0  TO  100 

P0VPYM=3  141593**2/ (  12.*EL0VD**2*EY) 

GO  TO  200 

100  POVPYMsSUOVS Y*  (  1  .  -  (EL0VD/K)**4/3.  > 

c 

200  WR I TE (  6  ,  lOOO  ) ELOVD  ,  POVPYM 


1OO0  FORMAT (//' *****»***FOR  L/D=  '.F6.3,'  MAX.  AXIAL  STRENGTH', 
*'  RATIO,  PU/PY=  '  ,  F  6  . 3 / / ) 

RETURN 

END 

C 

Cztxtxxztzzttttttzxzszxzxzzzxzxzxxxzzzzxztttztzxtxtztxtzzzxz 

C  .  .SUBROUTINE  TO  CALCULATE  MOMENT  FOR  ALL  SEGMENTS 

Czxxtxzszzxszxzzzzxzzxzstxtzxzzxxtzzxxxxtttxzzzxzzzzxzxzxzxt 

SUBROUTINE  MOMENT 


' 


' 


■ 


* 


■ 


. 


227 

228 

229 

230 

23  1 

232 

233 

23« 

235 

2  3  £ 

237 

2  3  £ 

239 

2  40 

24  1 

24  2 

242 

2  4  4 

245 

24  5 

247 

2  4  £ 

24  9 

250 

25  1 

252 

253 

254 

255 

256 

257 

2  5  £ 

25  9 

2  6  0 

26  1 

26  2 

26  3 

2  6  4 

26  5 

26  6 

26  7 

26  £ 

26  5 

270 

27  1 

272 

273 

274 

275 

276 

277 

276 

279 

280 

28  1 

282 

263 

2  84 

2  8  5 

286 

287 

2  8  £ 

2  8  5 

290 

29  i 

29  2 

293 

29  4 

295 

29  5 

297 

29  £ 

29  5 

300 

30  1 

302 

303 

304 

305 

306 

307 

308 

309 

3  1  0 

3  1  1 

3  1  2 

3  1  3 

3  1  4 

3  1  5 

3  1  £ 

3  1  7 

3  1  6 

3  1  9 

320 

32  1 

322 

323 

324 

325 

326 

327 

326 

329 

3  3  C 

33  1 

332 

333 

324 

335 

335 

33  7 

338 

339 
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I NTEGER  C , 0 . V 

COMMON  W|  lO)  ,W0  I  101  ,  EMOMYl  10)  ,  EMOMYTI  10)  ,  PDEL  I  1  O  )  .PHOPHYl  10)  , 

•SLOPE  (  1  O  >  ,  R  I  10)  ,  EY , ELOVD  ,  POVPY  ,  A  ,  PMUPH Y  ,  EOOQY  ,  EMOPY  (  101  , 

•EMOVP (  10)  , V , EMYOPY , NN , EMUT  , EMOMYM  <  4  0 )  .THETA0I4O)  .POVPYM,  JJ 
VsV  +  2 

EOOOYeFlOflTl V-21/IOOOO 
DO  1 OO  1=1,9 

EMOMYl  I  It6  *R  I  1  )*EOOOY*ELOVD 
EMOPY  (  I  )  =  EMOMYl  1  )  *  EMYOPY 
IF  (  EMOMY I  I  )  .  LT  .  O .  )  CO  TO  999 
100  CONTINUE 

RETURN 

9  9  9  WRITEIE  ,  1 OOO  » 

lOOO  F0RM4TI//1OX, 'NEGATIVE  MOMENT  VALUE') 

RETURN 

END 

C 

c 

Cssssssxs««x**ssstts(Stix»>*ts*t>*«tss*satzts*s>»<tSBSBSax(*t 

C.  SUBROUTINE  TO  CALCULATE  ESTIMATED  DEFLECTION  *  ELASTIC 

CsxatsaaaataxaasxataaaaaxaBsaataxaxataaaxtxsataaaxaaatattaaat 

SUBROUTINE  ELADEF 
I NTEGER  C , 0 , V 

COMMON  W(  10)  ,W0(  10)  ,  EMOMYl  10)  .EMOMYT  (  10)  ,  PDEL  I  lO)  .PHOPHYI  10)  , 

•SLOPE  <  10)  , R  I  10)  ,  EY  ,  ELOVD  ,  POVPY  ,  A ,  PHUPHY  ,  EOOOY  ,  EMOPY  (  10)  , 

•EMOVP  (  10)  ,  V  ,  EMYOPY,  NN,  EMUT,  EM OMYMl  40)  ,THETAO<40)  ,  POVPYM,  JJ 
C 

C.  CALCULATE  DEFLECTION  USING  EIGHT  DIVISIONS 
C 

THETA  =  SORT (  POVPY*  12  .  *  E  Y  )«ELOVD 

DO  1 OO  1=1,9 

EMOVP  (  I  IsEMOPYl  1  )  /  POVPY 

BETA  =  FLOAT  (  I  -  1  )*THETA/£ 

WO  (I  )  =  EM  0  V  P  I  I  )  *  <  (  (1  .  -COS  I  THETA  >  )  /  S  I  N  I  THETA  I  i*S]Nl  BETA  ' 

* ♦ COS  (  B  ETA  )  -  1  .  ) 

100  CONTINUE 
RETURN 
END 
C 

^sxxaaxxxxaaaxxxxfxxxtxxxxttxxxTxxxxxxxxxxxxxxaaaxxxtxxxxxxataxxxxxj 

C.  SUBROUTINE  TO  CALCULATE  CURVATURE  AND  DEFLECTION  FROM  M-P-0  RELATIONS 

[xxxxxxxxxxxxzxxxxxxxxxxxxxxxxtxxxxxxxxxxxtaxxxxxtxxxxxxxxxxxxxxxxxt 

SUBROUT  I  NE  PHIDEF 
INTEGER  C,C,V 

COMMON  W(  lO)  ,W0<  1 O  !  .EMOMYl  lO)  ,  EMOMYTI  10)  ,PDEL  (  1  O  )  .PHOPHYI  10)  . 

•SLOPE!  10)  ,  R  (  10)  ,EY,  ELOVD,  POVPY, A,  PHUPHY,  EOOOY  .  EMOPY'  10  , 

•EMOVP  (  10)  ,V,  EMYOPY,  NN,  EMUT,  EM0MYMI4O)  ,  THETAOl  40)  ,  POVPYM.  JJ 
C 

C.  .CALCULATE  VARIOUS  CURVATURES 
C 

DO  1 OO  N  N  =  1  ,9 
B  =  3  *  A  - EMOMYT INN  ' 

IF  IB. ECO.)  GO  TO  lOO 

PHOPH Y  I  NN  >  =  <  4  *  A  *  *  3  )  /  I  3  *  A - E M OM Y T I N N  )  ) *  *  2 

IFIPHOPHY(NN)  .GE  O  . AND  .  PHOPHY  (  NN  )  .  LE  .A'PMOPHYINN  I  a  EMOMYT  <  NN  1 
I F  I  PHOPHY  I  NN  )  .  GT  PHUPHY  ICO  TO  400 
1 OO  CONTINUE 
C 

C  .CALCULATE  SLOPES 

c 

SLOPE!  1  isPHOPHY  (  1  ) 

DO  200  J=  2  ,  5 
M=  J  -  1 

SLOPE! J)=SLOPElMi+PHOPHYiJI 
200  CONTINUE 

C  CALCULATE  DEFLECTION  ESTIMATES 
W  l  1  )  =0 
DC  300  K  =  2 , 9 
L  =K  -  1 

WIKieWlD+SLOPEILi 
300  CONTINUE 
RETURN 

400  1 F  I J J  .  GT  .  2  )  RETUR N 

WRITEIE  ,  2000 ) EMOMY  I  5  )  .EMUT 
JJsJJ-M 

2000  FORMAT  I  // 1  OX  ,  'MAX  CURVATURE  EXCEEDED  AT 

a'  M/MY=  '  ,  F  7  .  4  ,  BUT  ULT  .  MOMENT  =  ',F7  4) 

RETURN 

END 

C 

Caxxaaaaxxxtaxaaaiaxxxixaxaxtaxaxxxaaxxaxxx 

C  SUBROUTINE  TO  CORRECT  THE  DEFLECTIONS 

Qxxxxxaaxxxzxxxaaxzxxxxxxaxxxatxaxasaxaxxxx 

c 

SUBROUTINE  CORDEF 
I NTEGER  C , 0 . V 

C  OMMO  N  W I  1 O  )  , WO  I  10)  .EMOMYl  10)  .EMOMYTI  10)  ,  PDEL  I  10)  .PHOPHYI  10)  , 

•  SLOPE!  10)  , R  I  10)  ,EY,  ELOVD,  POVPY, A,  PHUPHY,  EOOOY  .EMOPY  I  10)  , 

•  EMOVP  f  10)  ,  V,  EMYOPY,  NN,  EM UT,  EMOMYM!  40)  ,  THETAO I  40  i  ,  POVPYM,  JJ 
ERR  T=W<  9  ) 

C 

C..  APPLY  CORRECTION  APPROPRIATELY 

C 

DO  100  1=1,9 

W I  I  >  =  -  <W(  I  1  -FLOAT  I  1-1  )/8.»ERRT) 

C 

C  CONVERT  TO  DEFLECTIONS  IN  TERMS  OF  L 

C 

Wtl  )«W<  I  > • 2  *EY*EL0VD/(64  .  ) 

1O0  CONTINUE 

RETURN 

END 

C 

Caaaazaattaatataataasxzx 

SUBROUTINE  STAR 

Caaaaxaaaxxaaxaaaaxaaaxxx 

C 


2000 


WR  I TE  <  6  ,  2000 ) 

FORMAT  I  //5X,  '••*»«•**«**•*•»*«•**«««•****»•********»******'  . 


■ 


■ 

. 


' 


/- 


■ 
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Appendix  A3 


Method  3 


1 

2 

3 

4 

5 

6 

7 

8 

9 

1  O 

1  1 

1  2 

1  3 

1  4 

1  5 

1  6 

1  7 

1  £ 

1  9 

20 

2  1 

22 

23 

2  4 

2  5 

2  6 

27 

2  8 

29 

30 

3  1 

32 

33 

3  4 

35 

3  6 

37 

3  £ 

39 

4  0 

4  1 

4  ; 

4  3 

a  c 

45 

46 

4  7 

4  £ 

4  S 

5  0 

5  1 

5  2 

5  3 

5  4 

5  5 

56 

57 

5  £ 

5  9 

6  0 

£  1 

6  2 

6  3 

6  4 

6  5 

6  e 

6  7 

£  £ 

£9 

70 

7  1 

72 

7  7 

7  4 

75 

76 

77 

76 

7  S 

80 

£  1 

S  2 

£3 

o  f 


c 

THIS  IS  THE  MOMENT  MAGNIFIER 

VERSION  OF 

THE 

c 

c 

I NTERRACT I  ON  D  I  AGRAM (METHOD 

3 

’  • 

c 

c 

IT  CALLS  PLOTIT  TO  PLOT  THE 

1 NTER ACT  I  ON 

D 1 AGRAM 

c 

DIMENSION  PO VP Y ( 40 ), EMOVMY 

(40) 

c 

R  E  A  L  *  8  TITLE  (  1  1  )/'L/D  =  0.05 

,  *  L/D=3 '  ,  ' 

L  /  D  = 

5  '  .  ' 

L/D  =  8  '  . 

’  L  /  D  =  1  0 

c 

*  '  L  /  D  =  1  2  '  ,  '  L/D= 1 5  '  .  '  L / D  =  1  8  ' 

' L /D  =  20 '  .  ' 

L  /  D  = 

23  '  , 

'  L/D  =  25 

’  / 

c 

READ  I  5  ,  500  )  (TITLES  (  I  >  .  1=1 

1  2  ) 

c 

READ!  5 , 500)  (TITLES  (  ]  )  .  1=1 

3 

,16) 

c 

P  E AD  ( 5 , 500 )  (TITLEStl  »  ,  1  =  1 

7 

,  20  ) 

c 

500  FORMAT ( 1 2 A4 ) 

c 

PRINT  TITLE  AND  READ  INPUTS 

WR  I TE ( 6  ,  2000  ) 

1  R EAD  '  5  ,  1 OOO ,  END  =  9 9  9  )  E L 0 VD , 

ALPHA . EUC 

M  =  0 

E YC= ALPHA* EUC 

c 

C.  .  PRINT  PARAMETERS  AND  TABLE  HEADS 
WR  I T  E  I  £  ,2100  t  ELOVD ,  ALPHA , EUC 
WR I TE  (  6  ,  2200  ) 

C 

C. COMPUTE  THE  DIFFERENT  LOADS 
IF  (  ELOVD  L  E  . O  .  )  G  0  TO  1 
P  1  =  3  1  4  1  S  9  3 

P I  SO: p I  *  *  2 
ND  =  1 
C 
c 

10  DO  1 OO  1=5,105,5 

POVPYIND  isFLOfiTI  1-51/100. 

C  CALCULATE  MO/MY 

POVPE  =  POVPY  (ND  >  *  12  *EUC*ALPHA*  <EL0VD**2  I/P1S0 
EMOVMY(ND)=( 1 . -POVPE ) * ( 1 . -POVPY (ND* ) 

C 

C  OUTPUTTHE  RESULTS 
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c 


c 

c 

c 

c 

c 


c 

c 


c 

c 


IFIEMOVMYIND •  IT  O  )G0  TO  150 
WRITE  I  6  ,  23OO1P0VPYIND  >  ,  EMOVMV  (  N  D  ' 

ND  =  ND  +  1 
1 OO  CONTINUE 

120  IF(POVPYIND) .LT.O  )G0  TO  200 
150  POVPY  I  ND  IsPOVPY ( ND  '  *0.001 
. . CALCULATE  MO/MY 

EMOVMY  I  ND  I  1  -POVPE  >*  I  1  .  -POVPY(ND)  ) 

IF  ( EMOVMY ( ND ». L T  O  > GO  TO  200 

.OUTPUT  THE  RESULTS 

160  WRITEI6  ,  2300  I  POVPY  I ND  )  ,  EMOVMY  (  ND  > 

200  CALL  CGPLIEMOVMY , POVPY , EMOVMV , ND , J , 1,1,3,0,0.0,0.25, 
*5.0,0  0,0  15,8.  .TITLES,  6  » 

200  CALL  PLOT  I T  <  EMOVMY  .  POVPY  ,  ND  .  J  ,1,  1  ,  3  ,  O  0,0.25, 

*5  ,0  ,0  15.6  ,6' 

CALL  CGPEP5(0.0,0  O.TITLEI J) ,8,0.1,0.75) 

JsJ+1 

.GO  TO  READ  FRESH  INPUTS 
GOTO  1 

998  WP  ITE  (  6  ,  2400  I  POVPY  (  ND  )  .EMOVMYIND  i 
9  9  S  J  =  O 

CALL  PLOTIT (EMOVMY , POVPY , ND , J , 1 , 1 ,3,0.0,0.25. 
*5.0,0.0,0.15,8  .6 

9  9  9  CALL  CGP L  ( EMOVMY  ,  POVPY  ,  EMOVMY  ,  ND  ,  0 ,  1  ,  1  ,  5 .  O ,  1  .  O  .  2  .  O  . 
*20  ,0,1.20,6) 

STOP 


C 

c  . 

C  FORMAT  STATEMENTS 
1 OOO  FORMAT (3F1C.5 ) 

2000  FORMAT  (  /// ,  23X  .  '  I NTER ACT  I  ON  CALCULATIONS  FOR  WOOD'// 
*23X,  'BE  AM  -COLUMNS-  RECTANGULAR.  SECTIONS'// 

*  3  4  X  ,  'WINTERS  (METHOD  3  l  4  /  / 

*23X,  '******************************»***'  ) 

2100  FORMAT  i  ///  10X  ,'  L /D=  *  .  F £  .  5 , 5 X ,  ' A L P H A =  '  .  F 8  5  .  / /  1 O X  .  ' E U C  = 

*  F  £  5  /  /  ) 

2200  FORMAT  t 10X, ' P/P Y ' , lOX, 'MO/MY'//) 


2300 

FORMAT  (  1  OX  ,  F5 

3  , 

1  2  X 

,  F  6 

4 , 1 1 X . F7 

4  /  /  ) 

2400 

FORMAT  (  / /  1  OX  , 

'  NEGAT 

I  VE 

VALUE  . . 

P/PY 

* 

'AND  M/MY  =  ' 

,  F  8 

4  / 

/  } 

STOP 

END 


file 


. 


. 
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Small  Scale  Specimens,  30 

Small-Scale  Tests,  43 

Strain  Distribution,  49 

Test  Procedure,  32 

TEST  RESULTS  AND  DISCUSSION,  43 

Test  Set-Up,  30 

Test  Specimens,  29 

Typical  Failure  Modes,  50 

Ultimate  Strength  of  Beam-Columns,  50 
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